Interlayer interaction and electronic screening in multilayer graphene 
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The unusual transport properties of graphene are the direct consequence of a pecuhar bandstruc- 
ture near the Dirac point. We determine the shape of the tt bands and their characteristic splitting, 
and the transition from a pure 2D to quasi-2D behavior for 1 to 4 layers of graphene by angle-resolved 
photoemission. By exploiting the sensitivity of the tt bands to the electronic potential, we derive 
the layer-dependent carrier concentration, screening length and strength of interlayer interaction 
by comparison with tight binding calculations, yielding a comprehensive description of multilayer 
graphene's electronic structure. 

PACS numbers: 



Much recent attention has been given to the electronic 
structure of multilayer films of graphene, the honeycomb 
carbon sheet which is the building block of graphite, car- 
bon nanotubes. Ceo, and other mcsoscopic forms of car- 
bon [l[ . Recent progress in synthesizing or isolating mul- 
tilayer graphene films has enabled access to their 
physical properties, and revealed many interesting trans- 
port phenomena, including an anomalous quantum Hall 
effect 
ature 



J , ballistic electron transport at room temper- 
7 1 , micron-scale coherence length 0, [1] and novel 
many-body couplings These effects originate from the 
effectively masslcss Dirac Fcrmion character of the carri- 
ers derived from graphene's valence bands, which exhibit 
a linear dispersion degenerate near the so-called Dirac 
point energy, Eu [lo| . 

These unconventional properties of graphene offer a 
new route to room temperature, molecular-scale electron- 
ics capable of quantum computing Q. For example, a 
possible switching function in bilayer graphene has been 
suggested by reversibly lifting the band degeneracy at 
the Fermi level (Ep) upon application of an electric field 
[ill , [l^ . This effect is due to a unique sensitivity of the 
bandstructure to the charge distribution brought about 
by the interplay between strong interlayer hopping and 
weak interlayer screening, neither of which arc currently 
well- understood. 

In order to evaluate the charge distribution, screen- 
ing, stacking order and interlayer coupling, we systemat- 
ically studied the evolution of the bandstructure of one to 
four layers of graphene using angle-resolved photoemis- 
sion spectroscopy (ARPES). We demonstrate experimen- 
tally that the interlayer interaction and stacking sequence 
affect the topology of the tt bands, the former inducing an 
electronic transition from pure-2D to quasi-2D when go- 
ing from one layer to multilayer graphene. The interlayer 
hopping integral, screening length, and layer-dependent 
carrier concentration are determined as a function of the 
number of graphene layers by exploiting the sensitivity 
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FIG. 1: (color online) Photoemission images revealing the 
bandstructure of (a) single and (b) bilayer graphene along 
high symmetry directions, F-K-M-F. The blue dashed lines 
are DFT bandstructure of free standing films [ij- Inset in 
(a) shows the 2D Brillouin zone of graphene. 



of TT states to the electronic potential. 

The films were synthesized on n-type (nitrogen, 1 x 
10^8 cm-3) 6H-SiC(0001) substrates (SiCrystal AG) that 
were etched in hydrogen at 1550 C. Annealing in a vac- 
uum first removes the resulting silicate adlayer and then 
causes the growth of the graphene layers between 1300 
to 1400Cjl3j- Beyond the first layer, the samples have a 
~ 1 monolayer thickness distribution; the bandstructures 
of different thicknesses were extracted using the method 
of Rcf. [llj. ARPES measurements were conducted at 
the Electronic Structure Factory endstation at beamline 
7.01 of the Advanced Light Source, equipped with a Sci- 
enta R4000 electron energy analyzer. The samples were 
cooled to ~ 30K by liquid He. The photon energy was 
94 eV with the overall energy resolution of ~30 meV for 
Fig.[2ra-d). 

The bandstructures of a single (Fig. [1] (a)) and a bi- 
layer (b) of graphene are reflected in their photoemission 
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FIG. 2: (color online) (a-d) The tt and tt* bands near Ef for 1 to 4 graphene layers, respectively. fc||= —1.703 A ^ corresponds 
to the K point, the corner of the hexagonal Brillouin zone. The F point is at fcy =0 A~^, while the M point is at -2.555 A~^. 
The dashed lines are from a calculated tight binding bandstructure, with band parameters adjusted to reproduce measured 
bands. Dark and light blue lines are for Bernal-type (ABAB and ABAC) stackings, while green lines are for rhombohedral-type 
stackings. (e-h) Photoemission intensity oscillation of tt bands at Ep — l eV as a function of fcy and k± momentum for 1-4 layers 
graphene. The photoemission intensity is normalized by the angle integrated intensity between Ep and Ef — 1.5 cV for each 
photon energy. 



intensity patterns as a function of fc|| . The data are com- 
pared with the scaled DFT bandstructures (see below) 



of free standing graphene layers IJ] . We identify one tt 
and three a bands, with the tt band becoming split in the 
bilayer. At the K point, the constant energy contour for 
the TT band becomes point-like near Ey P , 11 , 1^ , similar 
to graphite at the bulk H point [H, 1171 1^ Other features 
represent either underlying substrate states (especially 
in Fig. [Ija)) or replica bands due to diffraction from the 
underlying SiC surface reconstruction 0, [ll|, [l3| . 

The DFT bands are shifted to account for the Fermi 
level position, and expanded by 13% in energy to match 
the experimental total band width. This scaling effec- 
tively incorporates many-body interactions not included 
in the theory, as shown in an earlier ARPES study of 
graphite [l^. The experimental bandwidths of the films 
are in close agreement with those of graphite [l^ and 
thick graphene multilayers [ist . 

The TT bands near Ed exhibit a complex structure due 
to interlayer interactions ll| . Detailed photoemission in- 
tensity maps near Ep are shown in Fig. [^Ja-d) for N~l 
to 4 graphene layers. The measured tt bands are sup- 
pressed on one side of the Brillouin zone due to interfer- 
ence effects between the two equivalent sublattices [2o| . 
The number of tt bands increases with the number of 
layers due to interlayer splitting, clearly seen away from 
Ed , where the Coulomb potential of each layer does not 
play a major role (see below). The splitting between the 
highest and lowest tt bands increases with the number of 
layers and for quadlayer graphene (Fig [2l^d)) it is close 
to that of bulk graphite, ~0.7 eV [16|jl2l|. There is a gap 
between the tt and tt* bands in the bilayer (Fig. due 



to the inequivalent on-site Coulomb potentials in each 
layer [nj. 

The TT bands may be modeled by a tight binding (TB) 
calculation that takes into account the different stacking 
sequences and on-site potential ener gies , with a Hamil- 

13, 



Ionian generalized from Refs. 
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Here Ei is the on-site Coulomb energy for layer i, tt = 
Px + ipy, 7i is the interlayer hopping integral, v is the 
band velocity, and s=0 for Bcrnal (ABA. . . ) and 1 for 
rhombohedral (ABC. . . ) stacking. The Hamiltonian op- 
erates on the layer subspace i = (1,2,... N) while the 
2x2 operators a and P act on the {A, B) sublattice sites 
of the same or adjacent layers, respectively. The energy 
scale was defined such that E-p=Q and the Dirac crossing 
energy E'D^Tr TL/2N . The Hamiltonian can be readily 
generalized to arbitrary stacking orders (e.g. ABAC) by 
suitable rearrangements of the coupling operators Ps . 

It is well known that ABA and ABC stackings for 
graphene are energetically close, and stacking faults are 
commonly found in highly ordered pyrolytic graphite 
[2^. For samples with mixed stacking, we assumed 
the same on-site Coulomb potential for all stacking se- 
quences. For trilayer graphene (Fig. El^c)), we find two 
sets of TT bands resulting from different stacking se- 
quences. The blue and green dashed lines indicate TB 
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TABLE I: TB band parameters to reproduce measured band- 
structure for A'^=l-4 layers graphene (present work) and 
N=oo (graphite, Ref. [lal)- The electron density n is mea- 
sured in 10""^ electron per 2D unit cell. 
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FIG. 3: (color online) Potential and carrier concentration pro- 
files of the multilayer graphene as a function of the layer posi- 
tions. The electron potentials are shifted in the way that the 
potential of the outermost graphene layer is at zero. 



bands with Bcrnal and rhombohcdral stackings, respec- 
tively. The effect of different stackings is most apparent 
for the middle tt band near i^p— 0.5 eV, but is also seen 
in the upper and lower tt bands. 

The situation is somewhat different for the quadlayer 
shown in Fig.[2ld), where the measured and TB tt bands 
are compared. The dark and light blue dashed lines 
are for Bernal-type stackings, ABAB and ABAC, respec- 
tively. The green dashed lines are for rhombohedral-type 
stackings, ABCA and ABCB, which are equivalent within 
the nearest neighbor TB model employed here. We as- 
sumed that the on-site Coulomb potentials between lay- 
ers have the same sign and change monotonically. Four 
TT bands are well-reproduced by assuming Bernal-type 
stackings with the additional weak interspersed features 
suggesting minor contributions from rhombohedral-type 
stackings. The dominance of Bernal-type stackings for a 
quadlayer contrasts with the coexistence of Bernal and 
rhombohcdral stackings in the trilayer, and suggests the 
role of the second nearest neighbor in stabilizing Bernal 
stacking in bulk graphite. 

The tight binding parameters (Table I) reveal impor- 
tant details of the interlayer coupling and screening of 
graphene. Notice that our interlayer hopping integrals 
for multilayer graphene are significantly larger than that 
of bulk graphite, around 0.35 eV [16j, [21 1. In bilayer 
graphene, the interlayer hopping integral is reported to 
increase upon increasing carrier concentration We 
postulate that the larger interlayer hopping integral is 
similarly brought about by the higher carrier concentra- 
tions than that of graphite. 

The change in the out-of-plane coupling of the elec- 
tronic states is associated with the transition from pure- 
2D to quasi-2D as the number of layers increases towards 
eventual 3D (bulk) character. We study this transition by 
changing the electron wave vector normal to the graphene 
planes through a change of photon energy hv. Fig. [2je- 
h) shows the normalized photoemission intensities of tt 
bands at E-p — \ eV for 1-4 layers graphene as a function 
of electron momentum parallel (fc||) and normal (fcj.) to 
the surface, corresponding to hv ^ 17 - 210 cV. The 



maps in Fig. HJe-h) were normalized by dividing out the 
total intensity between E-p and E-p — 1.^ eV along M-K-F 
direction (Fig. [UJa-d)) at each hv. It is assumed that 
A:||= —1.703 A"-'^ at the K point with the same inner 
potential (~16.5 eV) as graphite 25, 26j . 

As expected for 2D systems, the tt bands display no dis- 
persion along A: (Fig.[2I^e-h)). For single layer graphene, 
the photoemission cross section decays smoothly and 
monotonically (this slow variation is normalized out in 
the data Fig. [2l^e)). For 2 to 4 layer graphene (Fig. [2jf- 
h)), however, the photoemission cross sections oscillate 
with a periodicity of ~2.0 A~^. This periodicity does not 
match with the reciprocal lattice vector of bulk graphite 
(0.932 A~^), which includes two graphene layers in the 
unit cell, but is close to the reciprocal interlayer distance 
of graphite (1.86 A^^). The slightly longer periodicity 
of the observed intensity oscillation may also imply a re- 
duced interlayer distance in our films, although precise 
assessment of the k± periodicity requires careful estima- 
tion of the inner potential. A reduction in the spacing 
for bilayer graphene was previously attributed to the in- 
creased screening with carrier concentration 1 1 1 1 . 

The intensity oscillations with k± for N> 2 are simi- 
lar to oscillations reported for quantized thin films states 
27, 28|, since the increase in the number of the features 



for 2 to 4 layers can be viewed as resulting from the 
quantum-well like nature of the wave functions in the 
graphene layers. The lack of such oscillations for mono- 
layer graphene results not just from the trivial lack of 
overlying graphene layers, but also due to a notable lack 
of coupling to the substrate, consistent with the lack of tt 
bands in the underlying SiC interface layer [l3|. This is 
in contrast to surface states lii-Eil, in which the photoe- 
mission cross section has a maximum near vertical tran- 
sitions of the bulk crystal, suggesting that single layer 
graphene wave functions have nearly ideal 2D character. 

The extracted on-site Coulomb potentials (Table I) 
may be used to estimate the screening length and the 
distribution of carrier concentrations. Fitting the on-site 
potentials (plotted in (Fig. OJa) with the outermost lay- 
ers' potentials aligned to the zero-reference level) to a 
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simple exponential decay allows direct determination of 
the interlayer screening lengths for three and four layer 
graphene, 1.45 A and 2.06 A, respectively. This trend of 
longer screening length for thicker films follows naturally 
from the decreased charge per layer, resulting in weaker 
interlayer screening. In turn, the estimated screening 
lengths are smaller than the reported value for graphite 
(3.8-5 A m,!!!), which has a much smaller carrier con- 
centration than our films. 

The total charge density n was determined by measur- 
ing TT-band Fermi surface areas, and was found to be al- 
most invariant with film thickness or stacking order (Ta- 
ble I); this explains the observed decrease of {E-p—Ed) 
with film thickness N (Fig. 2a-d). The ~10 % devia- 
tion of the carrier concentration for different thicknesses 
is due to the accuracy of the fitting procedure. 

Using Poisson's equation and the exponentially fitted 
profiles (Fig. [3lja)), the distribution of the measured to- 
tal charge density n (Table I) across the layers could be 
estimated for iV=3 and 4 graphene layers (Fig. [3I^b)). 
The carrier concentration decreases by about one order 
of magnitude for each adjacent layer. For the reported 
multilayer graphene devices the screening length is 
expected to be larger than the present case because of 
lower carrier concentration, therefore the carrier concen- 
tration profile is expected to be less steep. This suggests 
that the carriers in less-doped multilayers are distributed 
across several layers. 

While the electron potentials of a laterally confined 
surface state of vicinal Au(lll) were extracted using 
a combination of ARPES and scanning tunneling spec- 
troscopy 34|, |35[ , the evaluation of the potential, screen- 
ing length, and carrier concentrations solely from ARPES 
band structure measurements is unique to the present 
study. The present analysis is only possible because the 
topology of the tt states is very sensitive to the on-site 
Coulomb potential, and because of the very high energy 
and momentum resolutions of the experiment. 

In summary, we identify the tt band splitting due to in- 
terlayer interactions and different stacking sequences of 
the graphene layers. The interlayer interaction alters the 
character of the tt wave function from purc-2D in a sin- 
gle layer to quasi-2D in multilayer graphene. Exploiting 
the sensitivity of the tt states to the electron potential, 
the profiles of potential and carrier concentration in each 
layer and the screening length arc estimated for three 
and four layer graphene. We find the interlayer hopping 
integral and screening length deviate significantly from 
those of graphite because of the altered carrier concen- 
trations, which illustrate the unique electronic properties 
of graphene layers. 
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